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Cosmology with clustering anisotropics: disentangling dynamic and 
geometric distortions in galaxy redshift surveys 
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ABSTRACT 

We investigate the impact of different observational effects affecting a precise and accurate 
measurement of the growth rate of fluctuations from the anisotropy of clustering in galaxy 
redshift surveys. We focus here on redshift measurement errors, on the reconstruction of the 
underlying real-space clustering and, most importantly, on the apparent degeneracy existing 
with the geometrical distortions induced by the cosmology-dependent conversion of redshifts 
into distances. We use a suite of mock catalogues extracted from large N-body simulations, 
focusing on the analysis of intermediate, mildly non-linear scales (r < 50/;^'Mpc) and ap- 
ply the standard "dispersion model" to fit the anisotropy of the observed correlation function 
^(rx, rii ) . We first verify that redshift errors up to 5z ^ 0.2% (i.e. a, ^ 0.002 at z = 1) have 
a negligible impact on the precision with which the specific growth rate p can be measured. 
Larger redshift errors introduce a positive systematic error, which can be alleviated by adopt- 
ing a Gaussian distribution function of pairwise velocities. This is, in any case, smaller than 
the systematic error of up to 10% due to the limitations of the dispersion model, which is stud- 
ied in a separate paper We then show that 50% of the statistical error budget on p depends on 
the deprojection procedure through which the real-space correlation function, needed for the 
modelling process, is obtained. Finally, we demonstrate that the degeneracy with geometric 
distortions can in fact be circumvented. This is obtained through a modified version of the 
Alcock-Paczynski test in redshift-space, which successfully recovers the correct cosmology 
by searching for the solution that optimizes the description of dynamical redshift distortions. 
For a flat cosmology, we obtain largely independent, robust constraints on p and on the mass 
density parameter, Q.^ . In a volume of 2.4(/!^' Gpc)^, the correct Q.^ is obtained with ^ 12% 
error and negligible bias, once the real-space correlation function is properly reconstructed. 

Key words: cosmology: theory - cosmology: observations - large-scale structure of Universe 



1 INTRODUCTION 

The large scale structure of the Universe is one of the main ob- 
servational probes to discriminate among competing cosmological 
models and estimate their fundamental parameters, some related to 
space-time geometry, some related to the density fluctuations. The 
growth rate of density fluctuations, f( z), belongs to the s econd cat- 
egory. Since the pioneering works of lKaisen ( Il987h and lHamiltonl 
( Il998h . it was clear that one of the most promising ways to deter- 
mine f{z) is to exploit the apparent anisotropy in the clustering of 



galaxies induced by peculiar velocities, an effect commonly known 
as redshift-space distortions (RSD). 

Being f{z) directly sensitive to the mean density of mat- 
ter, for some time RSD have been used to est i mate the mass 
densit y parameter Hm (e.g. Peacock et alj 1200 Ic iHawkins et al 



2003: lda Angela et alj2005ijrbl:lRoss et al.l2007l : llvashchenko et al 
i2010i) . Later on, with the advent of other, more precise methods 
to estimate fljvi , like the Barionic Acoustic Oscillations (BAO), 
RSD began to be considered as a s ort of "noise" to be marginal- 
ized over (iSeo & EisensteinI I2003L 120071) . New interest on RSD 
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arose when it was realized tliat, if not marginalized over, they could 
tighten constraints over cosmological parameters (Amendola et al. 
2005). and in par ticular when it was shown (Guzzo et al. 2008 : 



Zhang et al.ll2008f) that they could represent a formidable tool to 



discriminate between a dark energy (DE) scenario and a modi- 
fied gravity theory for the origin of cosmic a cceleration. A num- 
ber o f forecast papers followed rapidly (e.g. iLinderiEOOSi; IWand 
120081 : ISong & Percival 2009). as well as ap pli cations to existin g 
and new datasets JCabre & Gaztaiiagall2009allbl:lBlake et alj|201 ih . 
As such, RSD are now regarded as one of the most promising tech- 
niques to extract precise estimates of /(z) from future redshift sur- 
veys. Besides their use as a probe to constrain alternative gravity 
theories, RSD can also be exploited in many different contexts. 
For instance, it has been recently demonstrated that RSD robustly 
const rain the mass of relic cosmological neutrinos JMarulli et alj 
I2OII ) and co uld be used to detect interactions in the dark sector 
JMaruUi et al.ll2012l) . Moreover, they can be of some help in astro- 
physical contexts, e.g. to investi gate the dynamical p roperties of the 
warm-hot intergalactic medium (lUrsino et al.ll201 ih . 

To effectively discriminate among competing cosmological 
models, however, one needs to measure the growth rate with per 
cent accuracy. This goal has prompted several works aimed at better 
identifying and characterizing the sources or uncertainty. One as- 
pect of the problem is how to optimally infer the growth rate of fluc- 
tuations from the measured RSD quantities. The anisotropic pattern 
of galaxy redshifts in 3D space allows one to estimate the distor- 
tion parameter |3(z) = f{z)/b{z). Then, to obtain f{z), one needs in 
principle an independent estimate of the galaxy bias parameter b{z) 
which is ill-constrained by theory and difficult to meas ure from the 
data. For this reason, it has b e en suggested (see e.^. Iwhite et alj 
I2OO9I : IPercival & Whit3 120091: ISong & PercivaJI |2009|) to express 
the constraints in terms of the observed product P{z)Gf (z), which 
in the linear bias hypothesis equals the theoretical combination 
f{z)Of;{z)- In these expressions ol" and Og are the rms values 
in spheres of 8 h^^Mpc of respectively the galaxy counts and the 
mass. In this way, the values of P(z)o|'' (z) obtained from the 
data at different redshifts can be self-consistently compared to the 
curves /(z)og(z) predicted by the theory (once they are neverthe- 
less normalized to a reference Og as provided, e.g., by CMB mea- 
surements). 

Another aspect of the problem is the impact of the observa- 
tional errors, related to redshift measurements and to the nature of 
the objects used to trace the mass inhomogeneities. This is com- 
monly addressed usi ng the Fisher information matrix (IFisheil 19351 : 
iTegmark et al.lll997D . which has become the standard method to 
forecast statistical errors in the cosmol ogical parameters expected 
from planned redshift surveys (see e.g. Sapone & Amendola 20071; 
Wang 2008: Under 200 8: White et al. 2009: Percival & WhitL 
20091 : IWang et alJl2010iriA cauaviva & Gawisei j|2010l: iFedeli et^ 
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Maierotto et al.ll2012n . One limitation in the use of the Fisher ma- 
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trix is that it relies on assumptions (e.g. the Gaussian nature of 
errors) that in some applications might not be fully justified. In 
addition, the Fisher matrix formalism cannot say anything about 
systematic errors that may dominate the error budget in future ex- 
periments aiming at per cent accuracy. Finally, the use of the Fisher 
matrix is justified on large scales, where linear theory applies, but 
may fail on smaller scales because of the improper treatment of de- 
viations from linearity in dynamics and bias. In all these cases one 
needs to use numerical simulations. 

Early attempts to test the accuracy of the linear model of 



RSD (Kaiseg [l987|) t hrough numerical simulations and de velop 
improved corrections (IS coccimarroll2004t [Tinker et al. 2006) have 
been recently extended, following the renew ed interest on this 



topic (ITaruva et al.l201ll : Jennings et alJ2011al lbl: IOkumura & Jind 
l201ll : lKwan et alj2012l : ISamushia et all201ll) . The work presented 
here is part of this ongoing effort to understand the limitations of 
RSD estimators and bring this technique at the level required by 
precision cosmology. 

Finally, a potentially relevant source of systematic error 
in the measurement of RSD is related to intrinsic uncertainty 
on the background cosmology, which is needed when convert- 
ing redshifts into distances. This introduces further geometric 
distortions (GD) in the observed clustering pattern, which in 
principle can be used to estimate the background cosmologi- 
cal parameters through the so-called Alcock-Paczynski (AP) test 
(Alcock & Paczynski 1979). Several methods and attempts to im- 
plement the AP tes t using different datasets have been pro- 
posed t Alcock & Paczvnskil I1979I: iPhillippj Il994l : iRvderJ [1993 : 
iRvden & Melottlll996l : iMarinoni & BuzzilbOlOl). including those 



that look for anisotropie s in 3D clusterin g (Ballinger et al. 199g: 
Matsubara & Sutd 1 19961 : IPopowski et "aLl 1998: Hui et al. 19991 : 
Matsubara 2000; McDonald 2003; Nusser 2005i: lBarkanall2006l: 
Kim & Croflii2007. : iPadmanabhan & Whitai2008h . In practice, GD 
have a smaller amplitude than and are degenerate with RSD 
(Simpson & Peacockl uOlO"). such that only with very high quality 
data one can hope to disentangle the two effects. Earl y applications 
of the AP test thus failed t o get significant results ( Outram et al 
20041 : iHovle et alj I2OO2I : Ida Angela et all |2005a"i^, 'Ross et al 
20071), but the situation is rapidly improving ( Chuang & Wang 



2OI2I : iBlake et alj l201lh . It is fairly clear that GD will play 
a fundamental role in the estim ate of cos mological parameters 
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from future surveys (ISeo & Eisen stein 2003 
'20 lO": 'Kazin et al.''201^; iTaruvaet alj|201ll : 
Hawken et al. 2012), but also that, if not properly accounted for, 
they represent one further source of systematic error in the mea- 
surement of /(z) from RSD. 

In this paper we explore how the precision and accuracy of 
/ depend on the limitations of i) using observed quantities (an- 
gles and redshift) to infer cosmological distances in redshift-space, 
and ii) using them to construct and model two-point statistics. 
We therefore first assess the impact on the recovered P of red- 
shift measurement errors and of the intrinsic uncertainty on the 
real-space two-point correlation function due to the deprojection 
procedure. We then look in detail into the effect of GD and in- 
troduce a simple technique to isolate RSD. We show that, under 
the hypothesis of a flat background, this allows us to simultane- 
ously measure p and Q.m- Finally, we relax the flat background 
assumption and investigate how GD can be used to constrain the 
cosmological parameters that enter the Hubble function. These 
tests focus on intermediate scales (r < 50/?^' Mpc), where non- 
linear effects cannot be neglected. These are the scales where RSD 
have the highest signal-to-n oise, also in last-generation surveys like 
Wiggle Z jfTake et al.l20lTI), VIMO S Public Extragalactic Redshift 
Survey dOuzzo et al.llin preparationl VIPERS) and SDSS-III BOSS 
dEisenstein et al J2OI 11) . 

In a parallel, complementary work (^ Bianchi et al.ll2012h . we 
use the same numerical methods adopted here to study the depen- 
dence of the uncertainty on the measured growth rate on typical 
survey parameters, as the volume and the density and bias of the 
adopted tracers. 

The paper is organised as follows. In ij2]we start describing 
the exploited set of N-body simulations and how to construct mock 
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halo catalogues. Then we give a general overview on RSD and GD. 
In ij3]we investigate the impact of redshift errors, deprojection ef- 
fects and geometric uncertainties on the estimate of the RSD pa- 
rameters. The method to disentangle GD and RSD is analysed in 
21 Finally, in ^ we draw our conclusions. 



2 THEORETICAL TOOLS 

In this section we review the theoretical tools we use to derive the 
random and systematic errors in the estimate of f{z) from RSD. 
Since we are not concerned on galaxy bias, we will focus on the 
distortion parameter (3. 



2.1 N-body simulations and mock halo catalogues 

Since we are interested in RSD on intermediate scales, our analysis 
will rely on mock catalogues extracted from N-body simulations 
that fully account for non-linear dynamics. The need for a suffi- 
ciently large number of independent mock catalogues, each one 
with a volume matching that of the typical ongoing redshift sur- 
veys, imposes to consider very large computational boxes. 

The 'Baryon A coustic Simulation a t the ICC', BASICC, meets 
these requirements JAngulo et al.ll2008h . This simulation follows 
the dynamical evolution of 1448 dark matter (DM) particles with 
massMpart = 5.49x IO'^/z^'Mq in a box of 1340 h^ ' Mpc usiri g 
a memory-efficient version of the GADGET-2 code ( ISpringell2005n . 
The cosmological model adopted is a ACDM universe with Q.m = 
0.25, fiA = 0.75, h = tfo/lOOkms-'Mpc-i = 0.73 CTr = 0.9. A 
detaile d description of this simulation can be found in lAngulo et alj 
( 120081) . DM haloes have been identified by linking more than 20 
particles with a standard friends-of-friends algorithm, so that the 
minimum halo mass isMnjjn = 20-Mpait — 1.1 • 10 A^'Mq. 

Since we are not interested in realistic models of galaxies, we 
simply identify mock galaxies in the simulation with the DM haloes 
with M > Mjnin. Their effective bias is 



»efflZj 



fM„^^^^n{M,z)b{M,z)dM 
"}-n{M,z)dM 



(1) 



where n{M,z) is the halo mass function and b{M,z) is their linear 
bias. Our goal is to assess the errors on (3, so we need a reference 
value to compare with. This is obtained by dividing the expected 
value f{z) = £lj£ (z) by the effective bi as obtained from Eg . Q. 
As we have seen in the companion paper JBianchi et al.ll20I2n . the 
effective bias of our mock D M haloes is in good agreement with the 
model of lTinker et alj OOIOO . that we will consider as our reference 
bias model. For comparison, we will also consider the bias function 
predicted bv lSheth et al.l ( 120011) . The ~ 10% discrepancy between 
these two models gives an idea of current theoretical uncertainties. 
Because of the limited mass resolution of the simulation, we 
are forced to ignore substructures within haloes. As a result, the 
largest haloes should represent cluster of galaxies collapsed into 
single objects rather than individual galaxies. This limitation has 
the effect of underestimating the contribution of small scale pair- 
wise velocity dispersion to RSD. We can think of this undesired 
effect as an attempt of mimicking observational constraints, like 
the fact that fiber collision imposes a limit to the number of spec- 
tra that can be taken in crowded areas, or catalogue-making pro- 
cedures like that of collapsing clusters or groups of galaxies into 
a single object. Further analyses of mock surveys from the Mil- 
lennium run demonstrate that the conclusions of our study do not 



depend on the missing sub- structure in the BASICC DM halo cata- 
logues JBianchi etal.l2012h . 

Moreover, since we want to mimic redshift surveys designed 
to explore the Universe at the epoch in which the accelerated ex- 
pansion has started, we will mainly focus on the simulation output 
at z = 1. Outputs at z = 0.5 and z = 2 have only been considered 
to check the robustness of our results. To construct the mock halo 
catalogues, we consider a local (z = 0) observer and place the cen- 
tre of the z = 1 snapshot at the coiTesponding comoving distance 
Dc{z =1), where 



Dc 



' dz' 



^'^-i:i 



(4)' 



(2) 



and the Hubble expansion rate is: 

H{z)=H^^[slM^.^+zf+Sl^{ 



1+z) 



+Ode exp 



('r 



1+W(z) 

1+z 



0.5 



(3) 



fl); = 1 — f^M — O.DE, w(z) is the DE equation of state, and the 
contribution of radiation is assumed negligible. 

In doing so, we neglect structure evolution within the box. Un- 
less otherwise specified, our error analysis will be performed using 
27 independent mock catalogues obtained by dividing the compu- 
tational box in 3-' subcubes, each one about twice the size of the 
VIPERS survey. The mean comoving source density in our mocks 
is0.003(/7/Mpc)^ 

To specify the distribution of DM haloes in redshift-space, we 
take the comoving coordinates of each object and compute its an- 
gular position and observed redshift: 



Zobs 



^Zc + ^(1- 



c 



(4) 



where Zc is the cosmological redshift due to the Hubble recession 
velocity at the comoving distance of the halo, vii is the line-of-sight 
component of its center of mass velocity, c is the speed of light 
and a^, is the random error in the measured redshift (expressed in 
km/s). 



2.2 Modelling Redshift Distortions 

RSD are induced by galaxy peculiar velocities. On large scales, 
where peculiar motions are coherent, RSD will be different than on 
small scales, where the velocity field is dominated by incoherent 
motions within virialized structures. An effective way to charac- 
terize these distortions is by mean of two-point statistics, like the 
power spectrum or the two-point correlation function. In this work 
we focus on the latter. We will refer to the redshift-space spatial co- 
ordinates using the vector s, whereas we will use r to indicate the 
real-space ones. 

The estimate of the spatial two-point correlation function ^(r) 
is based on counting galaxy pairs separated by a relative distance 
r. To characterize RSD it is convenient to decompose the distances 
in two components perpendicular and parallel to the line-of-sight, 
r = {r^,r\\). In absence of peculiar velocities, the iso-correlation 
contours of ^(rx,n|)are circles in the {r^,r\\) plane. RSD mod- 
ify these contours in a characteristic fashion: for large values of r^ 
coherent motions squash the contours along the peipendicular di- 
rection, whereas for small r^ incoherent motions elongate the con- 
tours along the parallel direction, generating the so-called "fingers 
of God" effect (Jackson 1972i) . 
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In the linear regime, the velocity field can be determined di- 
rectly from the density field and the RSD amplitude is proportional 
to (3. In this limit and in the distant observer approximation, the 
two-point correlation function in redshift-space can be written in 
the compact form; 

^(^±,^||)lm = ^o(.0^o(A')+^2(i-)-P2(A') + ^4(^)/'4W , (5) 
where fj = cos 6 = su/s is the cosine of the angle between the 
separation vector and the line of sight, s = i/i'^^ +s}, and Pi are 
the Legendre polynomia ls ( Kaisei 1987 ; Lilie & Efstathioij|l989l : 
iMcGillll 19901 : iHamiltonl ll992l : lFisheret alj|l994h . The multipoles 
of ^(^x I ■5|| ) can be written as follows: 



^0=1 



4p 

3 



7 



5 



m-. 



Ur)-i.{r)] 



^4 



8p2 

35 



^{r) + h,{r),-h,{r) 



(6) 



(7) 



(8) 



where ^(r) is the real-space undistorted correlation function, 
whereas the barred functions are: 



^{r) 



^{r 



1 [dr^irY- 
r^ Jo 

r^ Jo 



(9) 



(10) 



Eq. lO is a good description of the RSD only at very large scales, 
where non-linear effects can be neglected. 

A full empirical model, that can account for both linear 
and non-linear dynamics, is the so-called "dispersion mod el" 
JPeacock & Doddslll99d : Feeble j|l98d : lbavis & Peebleslll983h in 
which the redshift-space correlation function is expressed as a 
convolution of the linearly-distorted function with the distribution 
function of pairwise velocities /(v): 



^(.s 



±,s\\) 



:|_%V/(V)^ 



'±,■511 ■ 



v(l+z) 



(11) 



lin 



where the pairwise velocity v is expressed in physical coordinates. 
In this paper, we test two different forms for /(v), namely 



/exp(v) 



1 



Ol2\/2 



exp 



and 



/gauss (v) : 



1 



oizv": 



:exp 



V2\v\ 
C?12 

^12 



(12) 



(13) 



JDavis & Peebi5ll983l : lFisher et al.|[l994l : IZurek et alj|l994 . The 
quantity a^ is independent of pair separations and is generally in- 
terpreted as the dispersion in the pairwise random peculiar veloc- 
ities. Here we rather regard it as a free parameter that quantifies 
the cumulative effect of small scale random motions and statisti- 
cal errors on the measured redshifts 5z, and focus on the distortion 
parameter (3. 

This simple model for RSD depends on a few quantities: two 
free parameters, p and 012, a reference background cosmology to 
convert angle and redshifts into distances and the true two-point 
correlation function of haloes (galaxies) in real-space, t,{r), that can 



be either derived from theory or estimated from the galaxy redshift 
catalogue itself. A theoretical expression for the galaxy coiTelation 
function can be obtained, for example, from the observed galaxy 
luminosity function in the framework of the Halo Occupation Dis- 
tribution by adopting a theoretica l prescription for the halo two- 
point correlation function (see e.g. lYang et al.ll2003h . Alternatively, 
it is possible to estimate t,{r) from the measured t,{s^,s«) , with a 
two-step deprojection procedure. First, the observed redshift-space 
correlation function is projected along .sii : 

S(a)^E(..x)=2|J" ds\%{s^,s\). (14) 

Then, the real-space conelation funct ion c a n be estimated 
from the Abel integral ( .Davis & Peeblesi Il983t ISaunders et al.l 
[l992): 



1 r" 

''II ■''' 



dr' 



, dE{/^)/dr^ 



(15) 



In this paper we adopt this second approach and assess the im- 
pact of deprojection by comparing the results with the ideal case in 
which we use t,{r) measured directly from the mock catalogues. 

To eva luate the correlation function in the simulation we use 
the lLandv & Szalav ( 1993) estimator: 



Ur) 



HH(r)-2HR{r)+RR(r) 
RR{r) 



(16) 



where HH{r),HR{r) and RR{r) are the fraction of halo-halo, halo- 
random and random-random pairs, with spatial separation r, in the 
range [r — 8r/2, r + 8r/2] and 8r is the bin size. Since we are inter- 
ested in estimating P at intermediate scales, we evaluate the corre- 
lation function in bins of size 5r = 1 /i^' Mpc out to 50 /j^' Mpc , 
both in the parallel and perpendicular directions. We have checked 
that pushing our analysis out to r --^ 100 /i^' Mpc does not change 
significantly the results. Finally, since to perform the deprojec- 
tion procedure one needs to specify the behaviour of ^(r) at small 
scales, we have linearly extrapolated L,{r^,r« ) in the range rli"'" < 1 

h^ ' Mpc . We have verified that the results do not depend on the ex- 
trapolation scheme adopted. 



2.3 Modelling Geometric Distortions 

To convert observed redshifts and angular separations into relative 
distances one has to assume a background cosmological model that 
does not, in general, coincide with the true one. This mismatch in- 
duces asymmetries or anisotropics that, if spotted, can be used to 
constrain the background cosmological model itself. This cosmo- 
logical test, commonly known as AP test, can be performed using 
the observed two-point coiTelation function. 

In absence of peculiar velocities, isotropy in galaxy clustering 
guarantees that iso-coiTelation contours are circles in the {r^,r«) 
plane if pair separations are estimated assuming the correct geom- 
etry. In this sense, t,{r^,r« ) can be considered a standard circle in 
much the same way as BAO and Supemovae are considered stan- 
dard rulers and standard candles, respectively. The choice of an in- 
correct cosmology will distort these circles in a way that we can 
accurately predict. Indeed, the relation between comoving separa- 
tions in two different geometries reads: 



'■±1 



Pa,i(z). 

Da.2{z) 



r±i; 



fiiiz) 



r\\2- 



(17) 
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Figure 1. The impact of redshift errors on the measured real- and redshift-space (angle-averaged) correlation functions. The upper sections of each panel 
show the real-space true correlation functions ^(r) (black curves), the redshift-space correlation functions ^(s) (blue curves) and the real-space correlation 
functions obtained from deprojection, ^oir) (green curves). Each line shows the correlation measured in one of the 27 independent s ub-boxes (see Sec tion[2). 
The central and m iddle sections of each panel show the ratios ^(i)/^(r) and 5(.v)/^o(r), respectively. Linear theory predictions of [Tinker et al.l ( 12010.) and 
ISheth et al.l 1200 ID are indicated by black solid and red dotted lines, respectively. The grey horizontal bands represent 10% theoretical uncertainties. The panels 
refer to different redshift en'ors, as indicated by the labels. 



where the subscripts 1 and 2 refer to the two cosmological models 
and Df^ is the angular diameter distance: 



Da{z) 



//o(l + z)%/=I\ 



sm[^-Q.i^Dc{z] 



(18) 



These relations can be adopted to perform the AP test us- 
ing the two-point coiTelation function as follows. In the quest for 
P, the choice of a background cosmology is made twice: first to 
estimate t,{s^,su)from observed redshifts and angular positions 
and then to model t,{s^,s«) . We shall call assumed and test cos- 
mology the two cosmological models assumed to measure and 
model ^(ixi'^'ll) , respectively. The test is performed by computing 
^a('?xi'^H) f°'' ^ given assumed cosmology and comparing it with 
the model t,t{s'^ , in ) estimated for a test cosmology and rescaled to 
the assumed one (Eq. lll7t): 



^a(i±,i||)=^1 



PA.t 

Da,, 



H, 
H, 



(19) 



The correct values of P, £)a and H are found when t,t{s'^,si) = 
'%-i(s''^,sV), within the errors. Unfortunately, the small amplitude of 
GD with respect to RSD and inaccuracies in modelling the lat- 
ter make results obtained through this procedure not very robust. 
Moreover, this method depends explicitly on the bias model as- 
sumed to derive the DM correlation function from the observed 
galaxy positions. In this work, we have developed an alternative 
procedure, described in Section |4] that exploits both RSD and GD 
directly, and does not require modelling the shape of the DM cor- 
relation function and bias. 



3 MEASURING p FROM RSD 

3.1 The impact of redshift errors and deprojection 

Random errors in the measured spectroscopic redshifts contami- 
nate the clustering signal at all scales in a way similar to that of 
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Figure 2. Iso-coiTelation contours of ^{i^ , in ) measured in the 27 mock catalogues. Contours ai'e drawn in correspondence of the congelation levels ^(i± , s» 
{0.1, 0.2, 0.5, 1.5}. Different panels refer to different amplitudes of the redshift errors, as indicated in the labels. 



random peculiar motions at small scales. So far they have not been 
considered in the error budget, as local surveys typically had fairly 
precise redshift measurements errors (< 100 km s^'). A careful 
assessment of their impact is now in order, especially in view of 
future surveys using different observing techniques. For example, 
Euclid will be measuring redshifts from single emission lines in 
sUtless spectroscopic observations, with a requirement on the er- 
rors o f O;. < 0.001(1 -|-z) , which corresponds to 600 km s^' at 
z = 1 l lLaureiisetaUl2011[) . 

To assess the impact of redshift eiTors on the estimates 
of P, we have perturbed the redshifts of the mock galax- 
ies (Eq. lU) by adding a Gaussian noise of amplitude o,, = 
{0, 200, 500, 1000, 1250, 1500} km/s . These values cover a range 
extending out to errors that get close to (yet not as large as) those 
from photometric estimates (which in the best cases have <J,. ~ 
9000(1 -\-z) km/s). In what follows we express these errors as per 
cent uncertainties 5z[%] = {0, 0.07, 0.2, 0.3, 0.4, 0.5} (Eq. igj). 

To focus on the impact of redshift errors we do not consider 
GD, i.e. we assume the correct background cosmology. Moreover, 
we restrict our analysis to the case of Gaussian redshift errors. 
This is a quite accurate assumption for modelling the redshift er- 
rors o f typical spectroscopic galaxy surveys (see e.g. Lilly et alj 
12009!) . Furthermore, even in photometric galaxy surveys, it has 
been demonstrated that the adoption of a Gaussian distribution to 
simulate the impact of redshift er rors represents a reasonable ap- 
proximation (see e.g. lCunha et al.ll2009; Saglia et al. 2012) . We do 



not consider the impact of the so-called catastrophic errors, the ones 
caused by the misidentification of one or more spectral features, 
for the following reasons. In the assumption that the catastrophic 
outliers represent a perfect isotropic population (i.e. if such errors 
have a flat distribution), their effect is to reduce the amplitude of 
the correlation function at all scales. So they do not induce addi- 
tional distortions in galaxy clustering and do not bias the estimate 
of p. On the other hand, clustering distortions might be generated 
by systematic misidentification of spectral features. Their impact 
on P can only be estimated with mock galaxy catalogues, taking 
into account all observational effects that may vary case by case. 



3.1.1 Real- and redshift-space correlation functions 

In Fig. [T] we show the redshift-space two-point correlation func- 
tions, ^(.?), measured in the 27 mocks (blue curves in the upper 
part of the six panels) and compare them to the real-space ones, 
t,{r) (black curves). Redshift errors suppress the clustering ampli- 
tude on progressively large scales, as expected. A first quantitative 
assessment of this effect can be obtained from the ratio ^(.s)/^(r) 
which, in the linear limit, is simply related to P: 






1 



2p 
3 



(20) 



We plot this ratio in the middle part of each panel in Fig. \T\ (blue 
curves). When 8z = 0, this ratio is constant for r > 3 h^^ Mpc . 
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Figure 3. Large upper panels: best-fit parameters p and On, as a function of the redshift error 8z. Small bottom panels: systematic and statistical errors, as a 
function of bz- Specifically, the upper half windows show the best-fit values of p, averaged over the 27 mocks, and their scatter, o(p), represented by the error 
bars. The lower half win dow s show the random component of the relative error, o(p)/p[%]. Black solid and red dotted lines represent theoretical expectations 
from lTinker et al.N2010l) and lSheth et al.|j200ll) . respectively. The grey bands represent a theoretical uncertainty of 10%. The assumption on the adopted shape 
of f{v) and on ^(r) are labeled in the top-left part of the large upper panels. Top-left part of the figure: f(v) exponential form and true ^(r). Top-right: f{v) 
Gaussian form and true ^(r). Bottom-left part of the figure: f(v) exponential form and deprojected ^(r). Bottom-right: /(v) Gaussian form and deprojected 



On these scales, the value of P obtained from Eg. ( I20t is consistent 
with th eoretical expectations of lSheth et al.l ( 120011) and lXinker et alj 
( 120100 , represented by the red dotted and black solid horizontal 
lines, respectively. The horizontal grey band is plotted for refer- 
ence and represents a theoretical uncertainty of 10%. Increasing 
redshift errors to 8z ~ 0.2% has the effect of suppressing the clus- 
tering amplitude on ever larger scales and reduces the range useful 
to measure (3 but does not bias its estimate. For 5z > 0.3% the ra- 



tio t,{s)/t,{r) is biased high in the ever shrinking range of scales in 
which this ratio is constant, hence inducing a systematic error on P 
obtained from Eq. ( I20t . 

As a further step towards a realistic estimate of p, we also 
assess the impact of the deprojection procedure described in Sec- 
tion[3] The ragged green curves in the upper panels of Fig. [T]show 
the real-space two-point correlation function obtained from the de- 
projection procedure, t,£). The corresponding ratio ^{s)/b,£){r) is 
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Figure 4. Upper panel: eifect of GD on the halo correlation function in 
real- (middle set of black curves) and redshift-space (upper set of blue 
curves). The curves in each set have been obtained varying I2m in the range 
0.2 ^ I2m SJ 0.3 (from top to bottom). The lower set of grey curves shows 
the DM correlation function computed with CAMB (Lewis & Bridle 2002). 
The central green lines refer to the choice of the correct value J2m = 0.25. 
Lower panel: mean fractional error on the estimated ^ as a function of i^M in 
real- and redshift-space. 



shown in the bottom panels. Interestingly, the deprojected correla- 
tion function is in good agreement with the true ^(r) even for large 
values of 8z, indicating that the deprojection procedure does not in- 
troduce significant systematic errors. However, it increases random 
errors represented by the scatter among the ^(.s)/^D(r) curves. 



3.1.2 ^from thefiilIfitoft,{s^,su) 

A better estimate of (3 can be obtained by comparing the measured 
^{s± ,s«) with the model described in Section[3] The iso-correlation 
contours of ^(jx , •S|| ) calculated in the 27 mocks are shown in Fig.|2] 
for different values of 8z, indicated in the panels. Contours refer to 
the iso-correlation levels ^(i'xi'^'ii) = {0.1, 0.2, 0.5, 1.5}. The effect 
of increasing redshift errors can be clearly appreciated. The case of 
no errors (5z = 0) is characterized by the expected squashing of the 
contours at large separations induced by coherent motions whereas 
on small scales the fingers-of-God elongation is hardly visible. As 
already pointed out, this is due to the lack of substructures in the 
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Figure 5. Iso-correlation contours for 5(rx,r||) (left panel) and 
^(sx,.Sj|) (right panel). Contours are drawn at the iso-correlation lev- 
els {0.025, 0.05, 0.1, 0.2, 0.5} indicated by the labels. Different contours 
drawn at the same congelation level refer to the different values of 
SIm considered. The green contours refer to the correct cosmological 
model, ^M=0.25. 



DM haloes, such that the velocity field within virialized structures 
is poorly sampled. When redshift errors are turned on, fingers-of- 
God distortions appear and dominate the distortion pattern out to a 
scale that increases with 5z. 

Comparing the correlation function "observed" from the 
mocks in Fig. |2] with the model presented in Section |Z2l constrains 
the free parameters P and O12. This is done by minimizing the stan- 
dard X^ function: 






yu 



-4i 



% 



(21) 



where v;/ = ^SL,i-,%j) and yj^ = ^*^(.sx.;,.s||j;p,cyi2) are the 
measured and model coiTelation functions, respectively, and 5,y = 
8^(j I i,'S^ || j) is the statistical Poisson noise estimated following 
IMo et alj j 19921) . In each of the 27 mock catalogues, we fitted 
over the range 3 < rx,n| < 35 /i^'Mpc, with linear bins of 1 
hr^ Mpcboth in the parallel and perpendicular directions. As we 
have explicitly verified, the results presented in this paper do not 
depend on the particular form of Eq. 1 1211 and on the definition of 
clu stering uncertainties 5^(ixi'5||) (for a more detailed discussion 
see lBianchietalJl2012l) . 

The results are summarized in Fig.H] Let us focus on the upper 
left part of the figure. The points in the top panel represent the best- 
fit values of p and 012 obtained from each mock catalogue. The 
different symbols and colours indicate different redshift errors 5z, 
as specified in the labels. The best-fit p values should be compared 
with theoretical expectations using the iTinker et al.l (120101) model 
(black solid vertical line), which, as discussed previously, is a very 
good description of the intrinsic linear bias of our simulated haloes. 
The Sheth et alj 1 1200 11) model (red dotted vertical line) is shown 
for comparison. The vertical grey band shows the 10% uncertainty 
interval. These results have been obtained by comparing data with a 
model in which we have used an exponential form for /(v) and the 
true ^(r) of the DM haloes in the N-body simulation. Systematic 
and statistical errors and their dependence on 8z are quantified in 
the bottom panel. In the upper part we show the mean value of the 
best-fit P obtained by averaging over the 27 mock catalogues and 
its scatter, cj(P), represented by the error bars. In the bottom panel 
we show the rand om component of the relative error ct(P)/P [%]. 

As shown in lBianchi et alj <2012) . such systematic error de- 
pends on the minimum mass (i.e. the bias) of the haloes considered 
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Figure 7. Upper panel: bias of the DM haloes in the simulation, b = 
(^halo/^DM)"^. as a function of separation, estimated for different values 
of Q.M in the range 0.2 Sj ^m ^ 0.3 (blue curves). Red dots highlight the 
case of a correc t Q.m = 0.25. Error bars represent statistical uncertainties 
(Mo et al 1 9921) . Expec ted values from the models of Tink er et al.l feOlOl) 
and Sheth et al.l 120011) are represented by the horizontal lines. Bottom 
panel: mean slope of b(r) in the range 10 < r < 30 /i^' Mpc . 



Figure 6. Upper panel: true correlation function in real-space (green dots), 
deprojected correlation function with Q.m =0.25 (black curve) and depro- 
jected correlation function with other values of Q.m (red curves). Bottom 
panel: mean fractional error < 8^ > as a function of Q.m ■ 



and tends to decrease up to masses of 10 M0. In Fig. [5] we see, 
however, that redshift errors larger than 5z '-^ 0.3% produce an op- 
posite effect, which cancels and then overcomes the intrinsic neg- 
ative systematic bias on p. Interestingly, the rms error remains in- 
stead substantially constant, when the real-space correlation func- 
tion is well known (upper panel, --^ 5% for the volumes considered 
here). 

The upper right part of the same figure shows the results ob- 
tained when we model the velocity distribution function f(v) with a 
Gaussian function instead of an exponential one. The main effect is 
a very significant reduction of the systematic errors. This is due to 
the fact that redshift errors, modelled as Gaussian variables, can be 
regarded as a random velocity field with a distribution function that 
obeys a Gaussian statistics. What we learn here is that when red- 
shift errors dominate over the pairwise velocity dispersions, then 
/(v) is best modelled by a Gaussian function with dispersion com- 
parable to 6z. This is demonstrated by the fact that, in the plot, the 
best-fit values of O12 are comparable to the amplitude of the input 
redshift errors when 5z is large. 

The plots in the bottom part of Fig. |3] are analogous to those 
shown in the upper half except for the fact that, in this case, we are 
considering the more realistic scenario in which ^(r) is not known 
a priori but obtained from deprojection. Uncertainties in the de- 
projection procedure increase random errors by a factor of 2-3, de- 
pending on the amplitude of Sz. 



3.2 The impact of geometric distortions 

Before looking in more details into how GD arising from the choice 
of a wrong cosmological background can actually be exploited to 
our benefit, we would like first to understand how they impact the 
measurements of the growth rate from RSD. We first investigate 
how GD affect the estimate of the correlation function and galaxy 
bias. We then focus on the measurement of |3. Specifically, we as- 
sume a flat cosmology (so that f2y\ = 1— Hm ) and investigate the 
effect of choosing an incorrect value of Q.yi in the range [0.2,0.3], 
in steps of AHm = 0.01. All the other cosmological parameters are 
kept fixed to their true values. For this set of experiments we set red- 
shift errors 5z = 0. Since the amplitude of GD is smaller than that 
of RSD, to appreciate their impact we need to minimize sampling 
errors, i.e. trace velocities and density fluctuations with as many 
haloes as possible. Thus, in the following we shall use the whole 
simulation box with all its haloes, instead of the 27 subsamples. 



3.2.1 Impact on the measured correlation function 

Fig. |4] shows the effect of GD on the measured two-point halo 
correlation function in real- (middle set of black curves) and 
redshift-space (upper set of blue curves). The lower set of grey 
curves represents the correlation function of the DM obtained by 
Fourier transforming t he mat ter power spectrum computed with 
CAMB (Lewis & Bridl3 l2002l) . which exploits the HALOFIT rou- 
tine ( Smith et alj|2003 ). In each set, the central, green curve refers 
to the correct choice of background cosmology, Q.yi =0.25. The 
other curves refer to Hm values ranging from 0.2 (top) to 0.3 (bot- 
tom). The choice of the incoiTect cosmology also distorts the shape 
of the computational box. To account for this spurious effect, the 
random objects used to compute ^(r) have been generated within 
the same, distorted volume. 
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Figure 8. Plot of the ratio between the redshift- and real-space correla- 
tion functions, computed from the simulation box after converting redshifts 
into distances using 10 different cosmologies, con'esponding to 10 values 
of Siivi in the range [0.2,0.3] (blue curves). The green dots mark the cor- 
rect value, fo r ^m =0.25, wherea s the horizontal lin es show theoretical pre- 
dictions from lTinker et al.l bOlOl) (black solid) and lSheth et"ai] i200ll) (red 
dotted). The grey band represents a 10% uncertainty around the Tinker et 
al. value. 



GD enhance/dilute the correlation signal on all scales and thus 
modify the amplitude but not the shape of the correlation func- 
tion. The effect, quantified by the width of each set of curves, 
is very small. It can be better appreciated in the bottom panel in 
which we plot the mean fractional residual of ^, < 5^ > [%] = 
((^(nM)-^(nM = 0.25)/^(£^M=0.25)), where the average is 
over the interval 1 < r < 50 A^' Mpc . Since to first-order GD do 
not modify the shape of ^(r), the value of < 5^ > quantifies the 
amplitude of the spurious boost in the correlation signal induced 
by GD. In correspondence to the values Sly^ =0.2 and 0.3, already 
almost excluded by current observational constraints, the boost is 

Fig. |5] shows the effect of GD on ^(r^ , ''|| ) (left panel) and 
^(^x , ^11 ) (right panel). Contours are drawn at the correlation values 
{0.025, 0.05, 0. 1 , 0.2, 0.5}. The different curves at a given correla- 
tion level refer to different values of Sl^ . The green contours refer 
to the true geometry. 

Fig.[6]demonstrates that GD have little impact on the deprojec- 
tion procedure. The green dots show the true correlation function 
in real-space. The black curve shows the deprojected correlation 
function obtained assuming the correct value of Slyi . The other red 
curves refer to the other choices of fljyi . The effect is the same as in 
Fig.|4] an incorrect value for Q.yi boosts up or down the correlation 
amplitude on all scales by a factor < 10% (bottom panel), similarly 
to when ^(r) is measured directly. 



3.2.2 Impact on the measured galaxy bias 

An interesting aspect of GD is that the two-point correlation func- 
tion estimated assuming an incorrect value fljvi is different from 
the coiTect two-point correlation function measured in a universe 
with fljvi. This fact has some practical consequences, for example 
in the measurement of galaxy bias. Estimates of galaxy bias can be 
obtained from the ratio of the galaxy and mass two-point correla- 
tion functions. For example, the bias of the haloes can be estimated 
as bb!io{r) = (4alo('-,fiM)/5DM('-,i^M))°-^, where ^halo('',^M) is 
the real-space halo correlation function measured assuming some 
value of fljyi and ^om(''!^m) is the mass correlation function in the 
same cosmology. We have computed fohalo('') for the haloes in our 
mock catalogues. Results are shown in the upper panel of Fig.l?] in 
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Figure 9. Per cent mismatch between the value of (3 (upper panel) and O12 
(lower panel) computed for different values of ^m and the one obtained as- 
suming the correct SIm : 8X = \X{Q.u)-X{Q.m = 0.25)]IX{Q.m = 0.25), 
where X can be either P or 012. The curves show the result obtained assum- 
ing the tiue ^(r) (solid, blue) or the deprojected one (dashed, red). 



which we show fehalo('') obtained for different values of SIm in the 
range 0.2 ^ Hm ^0.3 (blue curves, from bottom to top). Red dots 
refer to the correct cosmology a nd error bars rep resent 1 a statisti- 



cal uncertainties computed as inlMo et al 



I992h Hori zontal lines 
(l2010h and ls"hethetalj 



show the model predictions of [Tinker et al 
(2001). 

These results show that GD affect both the amplitude of the 
estimated bias and its scale dependence. To estimate the effect, we 
fit each curve in the plot with a power law function b{r) =A + j-r. 
The spurious scale dependence is quantified by the slope j that we 
plot in the bottom panel as a function of Sl^i . Ideally, it would seem 
possible to estimate Q.m by requiring that b{r) remains flat on those 
scales where it should be constant. However, the smallness of the 
effect and the theoretical uncertainties on galaxy bias prevent this 
technique to be applied to real data. 



3.2.3 Impact on the measured vahie o/|3 

To assess the GD impact on |3, we have repeated the same analy- 
ses presented in Section [3] i.e. we have estimated |3 from the ratio 
^{s)/i,{r) and by fitting the full^(jxi'S||) • The blue curves in Fig. [8] 
show the ratio between the real- and redshift-space correlation 
functions for 10 different values of Sly[ in the range 0.2 ^ Hm ^0.3 
(from bottom to top). The green dots refer to the true cosmol- 
ogy case. Error bar s show the statistical errors computed accord- 
ing to the Mo et al. ( 1992) prescription. Reference values accord- 
ing to Tinker et al . (2010) and Sheth et al. (2001) are shown by the 
black solid and the red dotted lines, respectively, with a 10% the- 
oretical uncertainty indicated by the grey band. The scatter among 
the blue curves is significantly smaller than theoretical uncertain- 
ties, thus indicating that estimates of P from RSD in the range 
3 < r[/!^' Mpc] < 50 are robust to the choice of Slyi . 
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Figure 10. Contour plots of the DM halo correlation function, ^(r^.Tn) (blue curves), in real (top panels) and redshift space (bottom panels), measured 
for Q.M = {0.1,0.2,0.25,0.3,0.4}. The iso-correlation contours correspond to ^(rx,r||) = {0.04, 0.06, 0.1, 0.2, 0.5}. The red lines in both sets of panels 
correspond to the model ^(i'x,.S]|) . In the top panels, they simply give the dynamically-undistorted real-space ^(r) computed for the given cosmology and 
replicated over Jt/2; in the bottom ones, they also include the RSD dispersion model (Eq. lilt ) for the best-fit values of P and C12 derived for that specific 
Q.M ■ The green dotted curves show the geometrically undistoried i,(r^,r\\) measured at the true cosmology (Q.yi = 0.25), for comparison. 



The alternative way to estimate p from ^(^x ; ■S|| ) confirms this 
result. Fig. |9] quantifies the amplitude of the effect. In the upper 
panel we show the per cent difference between P computed for a 
given Sly[ , indicated on the X-axis, and the one obtained assuming 
the correct model SIm =0.25. The solid blue curve refers to the case 
in which we model ^(.sx,.S||) using the true t,{r), whereas the red 
dashed curve shows the case in which ^(r) has been obtained by 
deprojecting ^(.SX)'S||) ■ In both cases, the impact of GD is rather 
small, especially if compared to that of redshift eiTors. The corre- 
sponding error induced on P is less than 1% over the whole range 
of SIm analysed, when the true t,{r) is used in the model. Using the 
deprojected t,{r), the maximum error on p rises to --^ 4%. Finally, 
the lower panel shows the impact of GD on the other parameter 
of the fit, the pairwise dispersion CJ12. The error on this parameter 
turns out to be larger than the one on P, rising to ~ 20% for the 
extreme values of Q.f^ . 



4 DISENTANGLING DYNAMIC AND GEOMETRIC 
DISTORTIONS 

In this section, we investigate the possibility to perform the AP test 
using ^(ix,'5||) . The goal is to constrain both P and the cosmolog- 
ical parameters that enter Eq. (|3) exploiting anisotropics in galaxy 
clustering induced by RSD and GD. 



4.1 The method 

As we have anticipated in Section |23I a common approach is to 
use Eq. ( |19l l to model the two-point correlation function in a se- 
ries of test cosmological m odels and compare it with the measured 
one ( IBaUinger et alj 1996). Here we adopt an alternative proce- 
dure, which is a generalization of the iterative method introduced 



in lGuzzo et alj 1 120081) and that is found to be robust. The method 
consists in repeating the measurement of the correlation function 
in different test cosmologies, and then modelling only its RSD. Our 
working hypothesis is that, by construction, the agreement between 
model and data will be maximum when the test cosmology coin- 
cides with the true cosmology of the Universe, i.e. without GD. 
The steps of the procedure can be summarized as follows. 

(i) Choose a cosmological model to convert angular positions 
and redshifts into comoving coordinates. 

(ii) Measure ^(5x,i'||) . 

(iii) Estimate the real-space correlation function, ^(r), required 
to model dynamic distortions (e.g. through the deprojection tech- 
nique). 

(iv) Model only dynamic distortions (e.g. through Eq. Ill) , and 
derive the best-fit values of p and O12 that minimize the x^ function 
given by eg. 1211 

(v) Save this specific minimum value of the x^, that we shall call 
F{{DA,H}i). 

(vi) Go back to point (i) using a different test cosmology and 
estimate a new value for F. 

Once the whole set of {Dji^,H}i has been explored, the "best of the 
best" set of parameter values (P,cTi2,Dai^) will be then identi- 
fied by the minimum value of F{{Dfi^,H}i). The main differences 
between this procedure and the usual one are that i) the observed 
and model correlation functions assume the same test cosmological 
model, and ii) once Dp^ and H are fixed, one only needs to model 
RSD. In the case of a flat ACDM background, the success of this 
strategy is guaranteed by the small covariance between Q.f^ and P 
(see e.g. IRoss et alj|2007h . As a consequence, one can obtain an 
unbiased estimate of ^m even for an incorrect choice of p and CT12. 
One advantage of our procedure is that it does not require the 
modelling of the galaxy bias. Since tha galaxy correlation func- 
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Figure 11. The "pseudo-likelihood" function AF as a function of ^m . 
obtained with the method described in Section 14.21 i) In real-space fix- 
ing p = ai2 = (black line). H) In redshift-space, but fixing p and C5\2 
to their best-fit values as derived with the correct cosmology (red line), in) 
Same, but leaving also p and an as free parameters, either assuming perfect 
knowledge of the true real-space ^(r) (blue line) or iv) recovering it through 
deprojection (green line). 



tion can be obtained directly from the data through the deprojection 
technique, it is not necessary to model the shape of the DM corre- 
lation function (at point (iii)). The only assumption of the method 
is the intrinsic isotropy of the clustering. 

One disadvantage is the computational cost, since one has to 
estimate ^(.s^ , su ) for each cosmological model to test. However, 
the use of optimized linked-list-, Tree- and FFT-based algorithms 
allows ^(ix5'S||)to be computed sufficiently fast, as to efficiently 
explore the parameter space without resorting to supercomputing 
facifities. Alternatively, instead of directly measuring the correla- 
tion function at different test cosmologies, it is actually sufficient 
to measure ^(i-^ , su ) in a fiducial cosmology and rescale the result 
to a test cosmology, using Eq. l |19l l. A second disadvantage is re- 
lated to the estimate of the eiTors. The best-fit parameters are found 
by minimizing a function F that does not obey a x^ statistics. The 
reason is that the data themselves, which in this case coincide with 
the measured ^{s^ > •*|| ) ' depend on Dp^ and H. Therefore, since the 
values of F evaluated at different Df^ and H do not refer to the same 
dataset, the function F does not follow a x statistics. As a conse- 
quence, errors on D^ and H have to be evaluated in a different way, 
as we shall see below. 



4.2 Joint constraints on Q.y[ and P 

We start our analysis considering the case of a flat ACDM model in 
which Hm , the mass density parameter, fully characterizes the ex- 
pansion history and geometry of the Universe. Fig.llOlillustrates the 
result of applying this procedure to the full catalogue of DM haloes. 
The different panels show the iso-correlation contours of the two- 
point correlation function measured in real- (black curves in the top 



panels) and redshift-space (blue curves in the bottom panels). Con- 
tours are drawn at the correlation levels {0.04, 0.06, 0.1, 0.2, 0.5}. 
Different panels refer to the different values of Hm used to compute 
distances and estimate the correlation function, as indicated by the 
labels. The green dotted curves are drawn for reference and show 
the predictions for the true cosmological model Sly^ =0.25. As such, 
in the central panel they coincide with the black and blue curves. 
The red curves show the coiTesponding model for the two-point 
correlation function obtained using the best-fit values of |3 and O12 
estimated at each value of Slyi . 

In real-space (top panels) the model for l,{r^,r\\)i?, simply a 
replica over all angles of the real-space correlation function ^(r) 
(i.e. no RSD are present, coiTesponding to setting (3 = CT12 = in the 
dispersion model). This is shown here to evidence the interplay of 
the two effects. It could be seen as an idealized case in which we are 
able to perfectly correct for RSD, or can hypothetically reconstruct 
the real-space galaxy distribution. The iso-correlation contours are 
thus circles in the {r^,r\\) plane, when the correct cosmology is 
used. The effect of GD when varying the cosmology is then quan- 
tified by the mismatch between the green and the black contours. 
As evident, the best-fit value for Sly[ can be found by minimizing 
the difference between the red and the black curves, which is in 
practice the AP test. The best agreement is found for Q.^ =0.25, as 
expected, showing that this procedure is unbiased. 

Similar considerations can still be applied to the redshift-space 
case (bottom panels). For a given Slyi , the amplitude of the mis- 
match is similar to that found in real-space. This fact validates the 
hypothesis that GD and RSD are substantially independent. The 
difference between red and blue curves is still minimized for the 
correct reference value, Sly^ =0.25, showing that the result is unbi- 
ased also in redshift-space. 

Let us then quantify the ability of the proposed technique to 
jointly estimate Slyi and p. As we described, the best-fit values for 
(P,cJi2,fiM) are found at the minimum of the pseudo-x function, 
F . Note that, in this procedure, both the model and the measured 
^(^Xjiii) depend on Sly[ . The same happens with the errors, since 
the number of pairs in each bin is modified by the presence of GD. 
However, we have verified that this effect is small and can be ig- 
nored. More specifically, the shape of the function F and the po- 
sition of its minimum are very insensitive to 5^(nM)- In Fig. II II 
we plot AF = F{Sly^) — F^\j^, where Fj^in is the minimum value 
of F found during the exploration of the cosmological parameter 
grid {Sly[ j in our case). As in Section [3] we fit over the range 
3 < r^,r« < 35 h^^Mpc, with linear bins of 1 /j^'Mpcboth in 
the parallel and perpendicular directions. 

To evidence how the technique operates in detail, and under- 
stand its possible limitations, we proceed in increasing steps, as we 
did in Fig.[TO] We therefore first test the validity of the best-fitting 
procedure in real-space, an ideal case that would correspond to a 
perfect subtraction (or absence) of RSD. In this case, the value of 
F quantifies the mismatch between the black and the red contours 
shown in the top panels of Fig.[TO] The corresponding function AF 
is represented by the black curve in Fig. [TT] The minimum of the 
curve is found for Q.m = 0.25, again showing that the fitting proce- 
dure gives unbiased results. 

The red curve refers to redshift-space. It shows the values of 
AF obtained as a function of Q.^ ,-, after fixing P and <Ji2 to the 
best-fit values computed at Q.m = 0.25. In practice, this is again 
an idealized case in which the correct distortion pattern is known 
a priori and used as a reference against the one observed when as- 
suming different cosmologies. The RSD model is, in other words, 
inaccurate for all choices of Q.m but for Q.m = 0.25. Also in this 
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Figure 12. The 68 and 95 per cent pseudo-likelihood probability contours in the Q.m — ^a (l^ft panel), Q.u — m'o (central panel) and wq — Wa (right panel) 
planes. In each case, the other two parameters not shown are fixed to their true values. These constraints have been obtained in redshift-space, using the true 
^(r) and fixing P and 012 to their best-fit values as derived with the connect cosmology. The red squares mark the correct cosmological parameters of the 
simulation. The dotted black line in the left panel shows the flat background case: Q.^ = 1 — i^M- 



case, the minimum of AF is found for Sly[ = 0.25, thus indicat- 
ing that the switch-on of RSD does not bias our estimate. Interest- 
ingly, in this case the minimum is sharper than in the real-space 
case (black curve); this can be explained as due to the stronger 
constraints posed by the RSD pattern in the observed ^(rx,r||) , 
which reduces simmetry with respect to the simple real-space case. 
In other words, this is telling us that, if we were able to know red- 
shift distortions perfectly, e.g. from an independent measurement, 
then the constraints on the background cosmological parameters 
from an AP test would be more precise than those expected from 
the standard real-space geometric test. This is shown in this simpli- 
fied case by the fact that the red curve yields a smaller uncertainty 
on Slyi than the black one. 

In the general case, however, we do not know a priori the 
amount of redshift distortions and we would rather like to esti- 
mate also P (and CT12), together with Sly^ . The resulting constraint 
is shown by the blue and green curves. For the blue curve, we have 
assumed perfect knowledge of the real-space coiTelation function 
'%{r) (i.e. we have measured it directly from real-space positions in 
the simulation). We have already seen how crucial this is, as an in- 
gredient in the RSD model. Also in this case, the minimum is found 
at the expected value Hm = 0.25, although the fit is less constrain- 
ing because of the increased degrees of freedom in the model. The 
green curve, instead, depicts what happens with the same freedom 
in the model but in the most realistic case when one reconstructs 
t,{r) by deprojecting ^(.?x > ■*|| ) ■ Although the estimate of Sl^ is still 
almost unbiased (i.e. the systematic errors are smaller than the ran- 
dom ones), the minimum is now much shallower, thus indicating 
lower constraining power, i.e. a larger statistical uncertainty in the 
recovered value. These errors reflect the uncertainties in the depro- 
jection procedure, which are responsible for the scatter among the 
green curves in Fig.[T] 

As previously discussed, the function AF does not obey a x 
statistics and therefore we cannot use the curves in Fig. [TT]to de- 
fine confidence interval and estimate errors on Slyi . Ideally, one 
should repeat the analysis using many halo catalogues extracted 
from the simulation. However, in our analysis we have already con- 
sidered the whole computational box and would need to run more, 
independent N-body simulations, which are not available. We are 
therefore forced to evaluate errors using techniques that are typical 
of error estimates from observational samples. Specifically, we use 
the "block-wise bootstrap" technique: we divide the box into 27 



independent sub-boxes and build several boostrap samples, each 
containing 27 sub-samples selected at random, with replacement, 
from the original dataset. The l-o errors are then evaluated from 
the scatter on the rel evant quantities among the bootstrap samples 
JNorberg et ai]|2009l) . 

We apply this technique to quantify the uncertainties on our 
estimated values of Sl^ , P and 012 when using the procedure de- 
scribed in this section as applied in a realistic situation, i.e. redshift- 
space with free parameters (i.e. the cases of the blue and green 
curves in Fig.llOt. We obtain a value Q.m =0.24± 0.03, correspond- 
ing to a l-o uncertainty of 12%, when the true ^(r) is used in the 
RSD model (i.e. the black curve). When using the deprojected t,{r), 
the error on SIm grows up to --^ 40%. Still, no systematic bias is ap- 
parent. 

Finally, all these results have been obtained assuming no er- 
rors on measured redshifts. We checked directly that the impact 
of these errors on our conclusions is indeed negligible, as long as 
8z < 0.2%. However, when 5z > 0.3%, the resulting systematic er- 
rors on p do propagate to Q.m and can bias its estimate. 



4.3 Constraints on curvature and on the DE equation of state 

In this section, we investigate how GD can help in detecting pos- 
sible deviations from a flat ACDM scenario. Let us assume a more 
general DE model with equation of state: 

PDE Z 



WQ+Wa- , 

Pde 1 + z 



(22) 



JChevallier & PolarskillJOoIl: lLindeijl2003l) . In this case, the rele- 
vant cosmological parameters are: Q.yi , Q./^, wq and Wa (see Eqs.[3j- 

[HHH. 

As described in the previous sections, our AP test exploits 
only clustering distortions and does not consider the information 
encoded in the shape of the correlation function. The advantage is 
that our method does not depend on the galaxy bias model. The 
drawback is that with no constraints on the shape of ^(r), the above 
parameters are degenerate. This can be seen in Fig [12] that shows 
the 68 and 95 per cent pseudo-likelihood probability contours in the 
Q.M — f^A, ^M — WQ and wq — Wa planes. In each plot, the other two 
parameters that are not shown are fixed to the true values. The red 
squares mark the cosmological parameters of the simulation. These 
constraints have been obtained in redshift-space, using the true t,{r) 
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and fixing p and CT12 to their best-fit values as derived with the cor- 
rect cosmology, so that the red curve in Fig. Illlcorresponds to the 
pseudo-likelihood contours along the dotted black line in Fig [12] 
that illustrates the case of a flat universe (Q./^ = 1 — SIm)- The flat- 
ness constraint is almost perpendicular to the degeneracy in the 
Q.M — i^A> and is the reason that allowed to constrain Q.m in the 
previous section. 



5 DISCUSSION AND CONCLUSIONS 

In this work, we have investigated some relevant limitations exist- 
ing when using the anisotropy of galaxy clustering to measure the 
growth rate of density fluctuations, while accounting at the same 
time for the extra distortions induced by the cosmology-dependent 
mapping of redshifts into distances. More specifically, we have as- 
sessed the impact of different types of uncertainties, both observa- 
tional and theoretical, on the estimated values of P, the anisotropy 
parameter closely related to the growth rate. We have then tested 
how well, in presence of RSD, the correct underlying cosmology 
can be inferred. 

The main results of these analyses can be summarized as fol- 
lows. 

• The impact of Gaussian redshift errors on the estimate of RSD 
can be assimilated to a generalized small-scale Gaussian velocity 
dispersion, which can be quantified in terms of a single parameter 
analogous to the usual pairwise velocity dispersion O12. 

• In catalogues with volume, density and bias similar to the ones 
analysed in this work, we can estimate P from RSD with an accu- 
racy of 5 — 10%, regardless of the redshift errors. A general scaling 
formula for the statistical error on p as a function of the survey 
parameters is calibr ated and presented in the companion paper by 
lBianchieta"i]J2012h . 

• With typical spectroscopic redshift errors (a,. < 600 km/s), 
the anisotropy parameter P measured using galaxy-sized haloes is 
systema tically underes t imate d by ~ 10%. This is discussed in more 
detail in lBianchi et alj ( l2012n . where it is also shown that this sys- 
tematic error depends on the bias of the haloes considered. 

• Larger redshift errors (Oy > 1000 km/s ) introduce an opposite 
systematic bias in the estimate of P, if not modelled properly. This 
can be partly alleviated using a Gaussian model for the velocity 
distribution function /(v), rather than the exponential one. Note, 
however, that this may be influenced by the fact that a Gaussian 
distribution has been assumed for redshift errors (which is, in any 
case, a realistic choice for spectroscopic observations). 

• A key ingredient in modelling RSD in a sample is a good 
knowledge of the underlying real-space correlation function ^(r). 
Random errors on p are increased by a factor ~ 2 when b,{r) is 
obtained through the deprojection of the observed t,{s^,s«) , with 
respect to when using the correct t,{r). 

• GD arising from an incorrect choice of the background cos- 
mology affect both the measured correlation function and its model, 
and thus can impact the estimate of p. However, we have seen that 
this is very small, meaning that the value of p can be recovered with 
similar accuracy even assuming a wrong cosmological model. 

• GD have an impact on the estimated galaxy bias. The effect is 
to introduce a spurious scale dependence in the biasing function on 
those scales in which it is supposed to be constant. However, the ef- 
fect is very small and of the same order of theoretical uncertainties 
in current bias models. 

• We have implemented and tested an alternative procedure to 
perform the Alcock-Paczynski test from the observed ^(ixi^^ll) , 



measuring simultaneously P and the parameters that enter Eq. (|3). 
This is based on the (verified) assumptions that the effect of RSD 
dominates over GD and that the best match between RSD obser- 
vations and the RSD model is realized for the correct cosmology. 
We have shown that this procedure is robust and the results un- 
biased, in the case of a flat ACDM model. We give a first, ap- 
proximated estimate of the uncertainty that can be expected for 
O.M through a block-wise bootstrap resampling. In a volume V = 
2.4 • 10^(/i~' Mpc)^, we find that the expected errors on ii^ are of 
the order of ~ 12%, rising up to ~ 40% if the deprojected t,{r) is 
used instead of the true one. The results are very insensitive to the 
accuracy of the model used to describe RSD and to the magnitude 
of redshift measurement errors (up to 5z ~ 0.2%). Finally, we have 
investigated how GD can be exploited to constrain both the curva- 
ture of the Universe and the DE equation of state. 

In this paper we focused on the analysis of a simulation snap- 
shot centered at z = 1. Clearly, we could have analysed a corre- 
sponding box at z = 0, but preferred to focus on a redshift range 
which is becoming more and more important with ongoing deep 
surveys like VIPERS (which has an effective redshif t around 0.8 
and stretches out to z = 1.4 with its brightest galaxies dGuzzo et alj 
in preparation)), and with future larger surveys. Also, we concen- 
trated our analysis on intermediate scales, r < 50/i^' Mpc , where 
most of the RSD signal lies. These scales will remain important for 
these studies also in future surveys in which larger, even more linear 
scales will be surely better sampled, but nevertheless not sufficient 
alone for reaching the per cent precisions we are aiming for. 

Finally, we also limited our modelling to the simple disper- 
sion modeL_Weare_aware, as we show in our companion pa- 
per teianchi et alj|2012h . that this is not a fully appropriate de- 
scription of clustering and RSD on such mildly non-linear scales, 
when the precision on statistical errors becomes high. This is 
very probably at the origin of the observed ~ 10% systematic er- 
ror on the recovered P and significant work is b eing performed 
to improve it (see e.g. Ide la Torre & Guzzoll2012l and references 
therein). Still, in its simplicity the dispersion model performs sur- 
prisingly well when compar ed to much more complicated expres- 
sions (e.g. lSlake et al.ll201lh and delivers statistical errors compa- 
rable or smaller than those of more sophysticated non-linear cor- 
rections jde la Torre & Guzzoll2012ft . The fact that the impact of 
non-linear effects on estimated errors is quite limited is also sug- 
gested by the close similarity of the errors on P esti mated as in this 
paper to those predicted by a Fisher matrix analysis (IBianchi et alj 
[20l3). 
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